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ABSTRACT 


This thesis models induced magnetic fields from the motion of seawater in the earth’s 
magnetic field analytically and compares the results to arctic on-the-ice magnetic fluctuation 
measurements. The oceans have various types of internal motions, such as internal waves and 
turbulence. This motion of seawater, which is a conductor, in the earth’s magnetic field induces 
a current density. This current density, in turn, induces its own magnetic field. This thesis 
models internal waves and upper layer ocean turbulence analytically. The corresponding 
induced magnetic fields are calculated using a static form of Maxwell’s equations and parameters 
for the Arctic are inserted. Comparisons are made with measurements from the Arctic Internal 
Wave Experiment (1985). The predicted fields from internal waves have magnitudes that are 
measurable and of the same order of magnitude as ionospherically generated fields. The 
predicted fields from turbulence are several orders of magnitude smaller than ionospherically 
generated fields. Besides giving information about internal motions in the ocean, the seawater 


induced fields are a noise source in magnetic anomaly detection. 


lv 


ie 


ah? is ie 


Vie 


TABLE OF CONTENTS 


eee BUCH ON - Seok co Seek Aa ret) dete chisel eee ee 


Bae GROUP et, ch Gel det eek Se Ge Oe 


A. OCEAN ENVIRONMENT ...... .- 


Be INTERNAL WAVES ® e e e ® ® e e e e ® ® 


C. ARCTIC UNDER ICE SEAWATER TURBULENCE 


MAGNETIC SEGCNATURES @s. 22% © 6. 3 wo te se. % 


® ° 


A. PROPAGATION OF INDUCED ELECTROMAGNETIC FIELDS . 


B. INTERNAL WAVE SIGNATURE ...... . 
ties CAI @ Ne ey ree cs ln Oat. et Hox Ne, “ee 
Jn EMOUCeG Magnet 1G. Mend h: Ve. Ctra, is 

C. UNDER ICE TURBULENCE MAGNETIC SIGNATURE 
ieee GeOMet Inve, fo sc Sue Se i) ce ee. ee a ete 
2. Turbulence magnetic signature ... 


Doe MME RSURED. MACNEI EG SPECTRA eo Vo afar fo. 2 


15a chp 1 5) Gr ES Se em ae ae ae a Oe ee ee ee 


A. EVALUATION OF INTERNAL WAVE MAGNETIC 


MODEL e ® ® ® e e e ° ° e e e ° e e e ® 


FIELD 


® 


B. EVALUATION OF TURBULENCE MAGNETIC FIELD MODEL . 


Vee 


i> 


Li 


17 


Lay 


20 


34 


34 


34 


39 


42 


42 


48 


V. CONCLUSIONS 258. ce Fc 


VI. RECOMMENDATIONS... 6 cee 0 ci - 
VEIT. “APPENDE CES ice taco oop ioie ge es oe 
A. CALCULATION OF INTERNAL WAVE INDUCED MAGNETIC 


FIELD 2.0. 6 Se OU 6 
CALCULATION OF =TURBULENCE INDUCED MAGNETIC 
PED ee eee es 
INTERNAL WAVE VELOCITY POWER SPECTRAL DENSITY 
MEASUREMENTS 2950s Seuecmeees 6 oe et ee 
eB, TO Kenn ®, POWER SPECTRUM CONVERSION 


CALCULATION (oe ot een er cn es ee 


LIST ©F REFERENCES” .cuc <> crimes ett ten tee ce 


Initial Distribution est eee. en oe 


sal 


53 


35 


5) (i) 


oF 


60 


62 


67 


68 


71 


Table I: 


Table 


Es 


Wavelength Values ....... e 


Internal Wave Model Numerical Values 


(Continuum Method) OP ea Or ee ee ee ee ee 


Table 
(mode 
Table 


Table 


III: 


Internal Wave Model Numerical Values 


structure method) ...... «6 « « « 


IV: 


V3 


Signal Characterizations .... 


Turbulence Model Numerical Values 


7 doles b 


16 


45 


45 


46 


49 


Figure 1: Typical Ocean Thermocline, Halocline, and 
Pycnocline ~(Gross, 1972)". « . « 9s) ssc ™emet so oe 
Figure 2: Geometry of the Earth's Magnetic Field... 
Figure 3: Geometry for Calculation of Internal Wave 
Induced Magnetic Field = «2 «5. sss 9 cieceete ce 
Figure 4: AIWEX Horizontal Velocity Spectra (solid 
lines) ; Garrett-Munk Model Spectrum (dashed lines) for 
comparison. (Levine, et al., 1987) . ......-. 
Figure 5: Geometry for Calculation of Turbulence Induced 
Magnetic Field... 3 « “@- 8 66 ‘pepe ts ten uence 
Figure 6: AIWEX Power Spectra for magnetometer (Czipott, 
Poaney; 1985) 9 26 6 ws) 6) ese. 5 oon eels 
Figure 7: AIWEX Model Spectrum, CEAREX Model Spectrun, 
and Extrapolated Spectrum for AIWEX Magnetometer Data 
(Continuum Method) ee ee ee ee ee ee ee ee eek 
Figure 8: AIWEX Model Spectrum, CEAREX Model Spectrum, 
and Extrapolated Spectrum for AIWEX Magnetometer Data 
(Mode Structure Method) .........-. © « © « 
Figure 9: Turbulence Model and Extrapolated AIWEX 
Magnetometer Spectra ...... 6 © © « © © «© « « 
Figure 10: AIWEX Vertical Velocity Power Spectra 


(McPhee, 1989) ® e e ® e e ® ® ® e e ® ® is] ® ® 


sValslal 


18 


19 


24 


35 


41 


43 


44 


51 


52 


I. INTRODUCTION 

Seawater is a conductor situated in the earth's magnetic 
field. Motion of this conductor from surface waves, internal 
waves, turbulence, etc. produces a motional electromotive 
force. The resulting current density induces its own magnetic 
Greld. Understanding the production and propagation of such 
magnetic signatures provides information about the underlying 
seawater motion. Also, these induced magnetic fields are a 
noise source for other magnetic measurements, such as magnetic 
anomaly detection. 

This thesis constructed analytical models of these induced 
fields and compared the predicted values with experimental 
measurements. One model considered internal wave sources and 
another ocean upper boundary layer turbulence. The internal 
wave model had versions for the Arctic and lower latitude 
oceans. 

The internal wave model started with the mathematical 
derivation of a governing equation for seawater velocity from 
fluid mechanics. The model assumed no rotation, zero 
viscosity, incompressibility, and slow variation of mean 
density with depth. Separation of variables provided a linear 
differential equation for the amplitude of the vertical 
component of velocity. The corresponding solution and the 


separation of variables equations together provided an 


analytical model of seawater velocity, v. The velocity 
frequency dependence was then tailored to fit observed 
velocity spectra for either the Arctic Ocean or lower latitude 
oceans. 

Next, the induced magnetic field was calculated. The 
motion of the conducting seawater with velocity v in the 
earth's magnetic field B (assumed constant) induces an 
electromotive force and causes a current density J. This 
current density J induces its own magnetic field B'. Because 
the frequencies involved are on the order of 107° Hz or less, 
displacement currents were neglected and the inductive field 
B' calculated with the Biot-Savart Law. The total induced 
field B' at a field point was calculated assuming an internal 
wave of infinite y-extent with wavenumber in the x-z plane. 
Integration over all x and y, and over z within a depth range 
Peon, toms) sy Leldeds Bc The wavenumber dependence was 
removed by two different methods. In one, B' was converted to 
a power spectrum, transformed from a one-dimensional to a two- 
dimensional spectrum, and integrated over wavenumber from a 


minimum, k to infinity ctoegive 5 -«oe In the other, 


min! 
application of a boundary condition to the vertical component 
of velocity produced a dispersion relation between wo and k,. 
Choosing the lowest mode of oscillation, using the dispersion 
relation, and converting to a power spectrum yielded B'*(w). 


Data from the Arctic Internal Wave Experiment (AIWEX) in 


1985 provided a means to compare the model B'* power spectral 


density as a function of frequency with experimental results. 
The predicted magnitudes for the induced fields B'* were in 
the same range as extrapolated magnetometer measurements. The 
large spatial extent of internal waves produced signals of the 
same order of magnitude as ionospherically generated signals. 
The first k-dependence method yielded frequency responses that 
do not follow the de £2 dependence of the data. The second 
method does yielded the 1/f or 7 fe frequency dependence 
expected. The success of the second method implied that the 
modal structure of the internal wave field must always be 
utilized in calculating the induced magnetic fields. 

The turbulence induced magnetic field model had a similar 
formulation. It assumed the region of turbulent seawater 
resides in the upper boundary layer of the ocean between 
fepens Di wand D, “and to be three dimensionally isotropic 
within this region. The seawater velocity v was analytically 
modeled as a linear superposition of plane waves propagating 
in three dimensions with various wavenumbers and frequencies. 
This motion in the earth's magnetic field B resulted in the 
current density J, which was utilized in the Biot-Savart Law 
to calculate the induced magnetic field B' at a field point, 
again neglecting displacement current. The total field was 
found by integrating over all x and y, and over z between De 
and D. Keeping the principal term of the B' expression, the 
induced field components were integrated over wavenumber from 


a minimum, k EOeeinif i hac y. Then the horizontal field 


min? 


magnitude, B',, was converted to a power spectral density as 
a function of frequency. 

This power spectrum was compared to the AIWEX magnetometer 
data. The frequency dependence of the model agreed with the 
observed 1/f* dependence of the data. However, the predicted 
magnitudes were several orders of magnitude below the data. 


The AIWEX signal in the 10°'-10°% Hz range appeared to arise 


from fields generated in the ionosphere. 


II. BACKGROUND 


A. OCEAN ENVIRONMENT 

The conductivity of seawater principally derives from 
dissolved salts. In general, conductivity decreases with 
depth. An average value for ocean conductivity is 4 (S/m). 

Seawater density varies in the ocean. Temperature and 
salinity are the main governing factors. Density as a 
function of depth generally displays three zones (Gross, 
1971). In the surface zone, from the surface to 50-100 m, the 
density remains approximately constant with depth due to 
mixing of the upper layer of ocean by wave action. The next 
zone, from 50-100 m to around 500-1000 m, forms the pycnocline 
zone. Density increases rapidly with depth. Generally, the 
density follows a monotonic increase because the temperature 
declines with depth. The bottom zone, below the pycnocline, 
is called the deep zone. Here, density slowly increases with 
depth due to decreasing temperature. Fig. 1 (Gross, 1971) 
depicts these zones and their general temperature, salinity, 


and density variations. 


Thermocline Halocline Pycnocline 


_ Surface zone 
SE Pyenocline zone 
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= Deep zone 
ou 
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a _—_— 
Temperature (°C) Salinity Density 
(parts per thousand) (grams per cubic centimeter) 

Figure 1: Typical Ocean MThermocline, Halocline, and 


Pycnocline (Gross, 1972) 


The earth's permanent magnetic field permeates the oceans. 
Since the field is approximately that of a magnetic dipole, 
the direction and magnitude of the field vary across the 
oceans. It goes from basically horizontal with respect to the 
surface of the earth and magnitude of about 30,000 gamma 
(nanotesla) near the equator, to basically vertical with 
magnitude of about 60,000 gamma at the north and south 


magnetic poles. Seawater motions of the ocean within this 


field induce the electromotive force and the resulting current 


density J. 


B. INTERNAL WAVES 

A change in density with depth ( dp/dz < 0 ) allows 
buoyancy oscillations within the ocean known as internal 
waves. The momentum and continuity equations for fluids are 
used to derive an equation for the vertical velocity of these 
internal waves. 


Following Gill (1982), the fluid momentum equation is 


FO +2Qxu=-p"Vp-g + vV-u, (1) 


where u = fluid velocity = ui + vj + wk, O = system rotation 
angular velocity, and v = kinematic viscosity. Assume 
isentropic motion, i.e. no viscous effects (v=0) and no 


rotation (N=0). The mass equation of continuity is 


pre 4 Vu = 0. 


Assume density is a function of potential temperature, 6, and 
salinity, s, independent of pressure, p. Using the chain rule 
and taking ® and s constant with respect to time over a 


a@trrerential element 


Dp . 9p DO, dp Ds _, (3) 
De ob Dt Os DG 


The continuity equation then becomes the incompressible fluid 
COonadlt1on 


MG) = (Oe (4) 


Assume small perturbations in pressure, p', and density, p', 





= Dan De = + / ap, es (5) 
oO f p Pp’. Pp; dz 9P og: 
The fluid momentum equation becomes 
Ou = / / 
= -Vp! - p'gk. (6) 


jae 


Applying the convective derivative with average velocity v, 


the time derivative of density becomes 





Dp _ op Ps cee. (7) 
DE ay en ohe * "az ve 


The fluid velocity components u and v can be eliminated by 
first taking the partial derivative with respect to time of 
the continuity equation 


a ry pbc eh emailed ea 8 
Te as Gee Otoy "RISE on a 


and then using Equation 6 in component form to replace time 
derivatives of u,v with spatial derivatives of p'. The 


resents 


Ow _ @p’ | &p! 


° dzdt Ox? ay? — os 





Combining Equation [7] and the time derivative of the z- 


component of Equation [6] yields 


ew 2 eee 
ar2 pa ee pe lezoc « ioe 


where N is the Brunt-Vaisala frequency, given by 

















Noon eee, Sle (11) 
py, a2 
Taking the horizontal Laplacian of Equation [10] 
2 Zz Z h 
O O X ea Yw + ney , ES (oi jw Pee oe ne a ) Oo pt12) 
Gee “dx= “oy- Conair P, ex? dy? dzdt 


Substituting Equation [9] in the right hand side and 
rearranging yields an equation for w, the z-component of 


velocity, 


og 0° 3? 1 Oo 3 >, @ Ae - 
ap? ee) BN TCE oneal N nets eas 0. (13) 














Assume that the vertical velocity, w, varies with depth, 
Zz, much more rapidly than the density,p,. Then 


0 oa Ow 


1 a 
Bo 62 Gz az? il 





which is known as the Boussinesq approximation. 


Equation [13] then simplifies to 




















ae CC, OC, & 


GOt* sox" dy voz. Ox? ae i ye 


This equation describes the behavior of the vertical velocity 
component of an internal wave disturbance. 


Separation of variables can simplify Equation [15] 


W(X, V, 2 b= (Olay (eZ (2) Ge) (16) 


The X, Y, and T portions have the form 








2 
OX 4 kx = 0, (17) 
Cra. 
with the following solutions 
X = X,e7*"; y= ye; T =e oe (18) 
The z part becomes 
2 
OZ . x2(4-Nl2)*) 2-0, (19) 
Oz Ww? 


This equation possesses the Sturm-Liouville form. Thus, its 
solutions are orthogonal eigenfunctions 2Z(z) with real 
eigenvalues. The vertical velocity solutions then take the 


form 


10 


Wie. .G) =. 2) pa Keay eee (20) 


Since they are orthogonal, any internal disturbance can be 
expanded as a sum of these solutions. 


Bor -2(2Z) ) try .selutrons of the form 


Tileza=erz) | ere (21) 


Substituting this into Equation [(19] gives for w(z) 








d*w yee os NZ) 0) Dee re 
ee eee 3 eg may a w= 0. (22) 
The term with N(z)? makes this, in general, a nonlinear 


differential equation. 

To simplify this equation, approximate p(zZ) as a piecewise 
continuous function over adjacent regions such that dp/odz is 
a constant in each region. Then N(z) will be a constant in 
each region. Equation [22] can be solved for each region, and 
the solutions joined with boundary conditions specifying 
continuity of velocity. Alternately, approximate p(z2) by a 
superposition of linear slope segments each with dp/oz equal 
to a constant. Then N(z) 18S a constant for each segment and 
Equation [22] is again readily solved. These approximations 
are seldom applicable to the real ocean (Gill, 1982) for 


accurate velocity profiles. However, the resulting solutions 


get 


do have the general shape and the oscillating character of the 
profiles resulting from actual p(z) distributions. 


With N(z) equal to a constant, Equation [22] becomes a 





second order, constant coefficient, linear differential 
equation 
d*# aw , 
a— + pw = 0, 23 
dz? dz p se 
7 2 
where 2 Ky B = ee — Mie Re) - k2] ie 


The solutions of the characteristic equation are 


} 25 SSS See (NE = 07) (42 4 k2) (25) 


+ S w? A 


Define 6 as 


: (26) 


In a stratified ocean, the maximum oscillation frequency 
equals the buoyancy frequency N. Thus, wo < N , and 6 is real. 


For 0 < N, the general solution of equation [23] is 


R(z) = eM? 1c, ete + Cc e|-f8z] | (27) 
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Initial and boundary conditions determine the constants Cl and 


CZ. 
The vertical velocity component is then 
] + + - . (28) 
ie yee (7) en x A ae Oa yet ke wt) 
The incompressible continuity condition, Equation [6], 
determines the horizontal components, wu and v, 
7 k : 
Gy ane) ist Se) Woes) hee Le (29) 
Ka Ke Pk 
“, k 
VAX Meh Z7 b= ist - k, W] —— eilk:r- ot). (30) 
4 Kacehe, 


These analytical expressions for the seawater velocity now 


facilitate calculating the induced magnetic fields. 


C. ARCTIC UNDER ICE SEAWATER TURBULENCE 
In the arctic, turbulent seawater flow exists from below 
the ice (effectively zero depth) down to some depth, D, at 
which the motion transitions to internal waves. 40m is a 
typical value for D (private communication with T. Stanton). 
Utilize superposition of three-dimensional plane waves of 
varying wavenumber and angular frequency, w, to model this 


turbulent region 


13 


v= [v,(k,,o) "Lf + vi(k,,w)euwF + v.(k,,0)e™?F) e-tet, (31) 


Assume three-dimensional isotropy of turbulence to simplify 


the analysis 


v,(k,,@) = v,(k,,@) = v,(k,,0) = V(K,®). (32) 


Empirically model the velocity amplitude function v(k,) as 


i — —— (33) 


The power spectrum for this function is then 


ve 


View) 2a eee (34) 
(k2 + k?) (w%, + w?) 


The constant factors k, and ©, in the denomimator Keeper 


power finite as k and w approach zero. For large k, this 
exhibits the same k™~@ dependence as for internal waves 
(Garret, Munk, 1972). For large o,this exhibits thomman 


dependence characteristic of under ice internal wave fields 


(ONR Report, 1991). The velocity expression is 


CK cD) ee ee Oe (35) 
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III. MAGNETIC SIGNATURES 


A. PROPAGATION OF INDUCED ELECTROMAGNETIC FIELDS 

This model neglected attenuation due to the propagation of 
the electromagnetic waves through seawater, a conductor, 
because of the low frequencies involved and the corresponding 


long wavelengths. The electromagnetic waves are characterized 


by their wavelength, A. For a "good Peonductor’”,. 2 ike 
seawater, and a "good insulator", like air, the wavelengths 
are 
27 Ee 
LN eae FS) eet Aenaaieroe xe a (36) 
frp of 
where 
pu =4nx10-7 0=CoOnd= 42, c=3.0x10°=. (37) 
Ae m S 


Table I lists wavelength values for typical frequencies of 


Paeerest in internal ocean seawater motions. 


iB) 


Table I: Wavelength Values 


Freq fH2) A doaater (m) Xair (m) 8 
10 1570 3x10 
en 5030 cpl 
10° 15,700 3x10! 
10° 50,300 3x10!! 
oes 157,000 3x10! 


The induced magnetic fields will "propagate" with 


characteristic wavelength A up through the ocean and then 


seawater 


with the much, much longer wavelength dae above the surface. 

Characteristic ocean buoyancy frequencies are around 10% - 
10°* Hz. Since internal wave frequencies are always less than 
the buoyancy frequency, and the depth of the Arctic Ocean is 
around 3000 m (Dietrich, et al., 1980), the corresponding 
wavelengths of the induced electromagnetic radiation are much 
larger than the depth. Within the ocean and in the air above 
the ocean surface, conduction current, J, and not displacement 
current, e€dE/dt, creates the induced magnetic field. The 
field point is in the "near field" and propagation effects are 
negligible. Relevant turbulence frequencies range around 10°!'- 
10°37 Hz. We confined our interest to turbulence in the upper 
boundary layer of the Arctic Ocean. Since a nominal depth for 
this layer is 40 m, the electromagnetic radiation wavelengths 
involved are again much larger than the depth and fields above 
the surface are "near field" with negligible propagation loss. 
Thus, we neglected attenuation due to propagation in the 
seawater conducting mediun. For these same reasons, the 


sea/ice/air interfaces do not attenuate the induced fields. 
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The earth's magnetic field inside the ocean will possess 
the same value as in the air above the ocean. This follows 
from the assumption that the earth's field is a static (DC) 
field, and that the magnetic susceptibilities of the air and 
ocean are the same (i.e. both are non-magnetic). The 
equivalence derives from the resulting boundary conditions 


between the two media 


B = B 


air,tangential’ 


B = B 


sea, tangential sea, normal alr,normal' 


B. INTERNAL WAVE SIGNATURE 

1. Geometry 

The method of Crews and Futterman (1962) allowed us to 
calculate magnetic signatures. The current density followed 
from the seawater velocity in the earth's magnetic field. 
This current density was then utilized in the Biot-Savart 
method of calculating magnetic induction at a field point. 

Fig. 2 shows the geometry of the earth's magnetic field 
within the ocean (northern hemisphere). To simplify the 
calculations, choose the wave vector in the x-z plane. Then 
the y-components of the wave vector and the velocity are zero. 
Assume the wave extends infinitely in the ty directions. The 
seawater in the volume element dxdydz at the point (x,y,zZ) 
with velocity v(x,y,z,t) forms a source element for an induced 


magnetic field at the field point located at (0,y,-h). Fig. 


Jey; 


3 shows this geometry. "z" represents the depth below the sea 
surface. "h" represents the height of the field point above 


the sea surface. 


Figure 2: Geometry of the Earth's Magnetic Field 
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Figure 3: Geometry for Calculation of Internal Wave Induced 
Magnetic Field 
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2. Induced Magnetic Field 

Seawater behaves as a relatively good conductor, with 
conductivity of approximately 4 S/m. The motion of this 
conductor in the earth's magnetic field creates an 
electromotive force. This in turn leads to a current density 


given by 


J=0 (VXB) =o ( v,B, - v,B,) f. (39) 


Since velocity is confined to the x-z plane, the current 
density J is in the y-direction. 

This current density creates a secondary induced magnetic 
field. The Biot-Savart law gives the differential induced 


magnetic field due to a current element di 





dp = Hb, di XR 
2m R2 


' di = J dA. (40) 

The prime indicates the secondary induced field. Integration 
of Equation [40] gives the components of the total induced 
field. The assumption of infinite extent of the wave in the 
+y directions (i.e. symmetry about y=0) results in the y- 


component summing to zero. The x and z components are 


Bl = [[ 4B: = [Le  (h+2), (41) 


fe . “a 
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Substituting in the expressions for di, v, 3B, and R and 


integrating {see Appendix A] gives 


Pee BOC ape te 6 e (estible a (ke ibe 3. eo (kati8)2 eke ible , (43) 
ae ——_—__— + ——-.—_) cosycos$ + —sind (——___—- - ——.— 
Ps 2 = ! ( -k,+i6 i k,+16 ” > Re na -k,+16 k,+18 ) D, 


, o) BoBC eer er ectiile eke ible 38 ; eg okrtibe : e (ke ibe A (44) 
pee a amar neg ee ee aes K-16) Pe 
Here, 6 reduces to 
oes Z 
Saute eee (45) 
* 2 
W 
These equations are fairly complicated. However, they 


simplify somewhat when evaluated specifically for conditions 
im. the Arctic 


(46) 


Bae kK: ® = 90°: cosd = 0 sing = 1; D, = 0 


Re 


The induced magnetic field components become 











PoOBC pn. N2-w2 1 {-k.#i8)D 1 (-k,-481D 2w? 
i) eee Nese dwt eel « a pee ee x Sener. ene : 7 
- oo ae Dy : [4 Serres lame Knt-207) | (47) 
i Use a N2-2 1 (-k,9i8)D 1 (-k,-48)D 20? 
Ba es pg Pe iat | se! x - P ees a (48 
z ee ze ( ¢ Preis © ‘Sere K,(N?-207) » ( ) 


iBi=BiL= (p BoC) oie a e [ k, - e*"(-k,cos8D - 3sindD) }. (49) 


k2432, Kk, 


Inserting the expression for 6, these become 


oe 


ipl=Bley BoC oe kltie-tet w/N?-w? (1 - e* (cos ee OT ED z Nt=OF ee Nt-0" kD a ( ) 
k, N? w? w? oe 


Note that the x and z components are exactly the same except 
for the phase factor of i 


d 


(51) 


tla 


Bo=i Bi, with ise 


Thus, the x and z components are 90° out of phase, and the 
induced magnetic field vector rotates in the x-z plane. It 


has magnitude 


IB’ = |B] = [B2| (52) 


Important characteristics of B' are (1) exponential decrease 
in magnitude with height, h, above the internal wave field, 
and (2) increase in magnitude with the vertical height D of 
the internal wave field. 

These expressions contain angular frequency wo, wavenumber 
ky and) tine” te Comparison with experimental data often 
requires the total field power spectrum as a function of 
frequency. Thus, all of the wavenumber and time dependencies 
must be removed. 

The frequency and wavenumber dependencies of the velocity 
were empirically modeled based on actual measurements of 
velocity power spectra. The amplitude of the internal wave 


velocity, C, which so far has only been a multiplicative 


component in the equations, became a function of frequency and 


Pa 


wavenumber. This, in turn, added additional frequency and 
wavenumber dependence to the induced magnetic field B'. 

The Arctic Internal Wave Experiment (AIWEX), conducted in 
1985 in the Canada Basin of the Arctic Ocean (depth > 3500 m), 
measured horizontal velocity with an electromagnetic current 
meter at a depth of 100 m. Fig. 4 shows one such spectrum. 
Sampling rate limited the highest frequencies measured to 


around 6 cycles per hour (loco. HZ, eneuocal. Sbuovancy 


frequency. The spectra indicate a 1/f frequency dependence 
over the range 107' - 10 cycles per hour (29esclO > = 2 esa 
Hz) (Levine, et al.,1987). This frequency dependence was 


empirically modeled by 


es (“A" for AIWEX) . (53) 


For W >> W,, it behaves as 1/w. For wo < Wo its value stays 
approximately constant and limits the total energy in the 
spectrum as w approaches zero. 

The Coordinated Eastern Arctic Experiment (CEAREX) in 1989 
measured velocity power spectra at depths of 50, 100,150,200, 
and 250 m below the surface of the ice. The 50 m measurements 
were made with an acoustic doppler profiler and the other 
measurements were made with electromagnetic current meters. 
Appendix C contains reproductions of the resulting velocity 
power spectra. All plots showed an approximately Ie) £6 


frequency dependence in the range 0.1-10 cycles per hour 
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Figure 4: AIWEX Horizontal Velocity Spectra (solid lines); 
Garrett-Munk Model Spectrum (dashed lines) for comparison. 
(Levine, et al., 1987) 
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(2278x100 BOA Tomo Hz). Often this dependence flattened 
out for frequencies below about 4x10° cycles per hour 
(1.11x10°° Hz). (Czipott and Podney, 1989). 


This frequency dependence was empirically modeled by 


‘ap ‘a 
lv.|? = ———- = ——_*+—_., (“C" for CEAREX) . 
Grea Ge (0 oma es ey 


(54) 


Bom .>> i > At behaves as Wo ANGE tO at <6, it remains 
approximately constant. Choosing , corresponding to 4x10 


cycles per hour gives 


@, = auf, = 2m (1.11x10”" Hz) = 7.0x10"' rad/sec. (S55) 


Measurements of internal wave spectra in other oceans of 
the world at lower latitudes have indicated a characteristic 
Wie frequency dependence. Garrett and Munk (1971) formulated 
a popular empirical model based on these observations with 
Such a frequency dependence. However, the AIWEX spectra 
revealed less total energy in the arctic internal wave field 
by factors of 15-30 than in fields described by the Garrett- 
Munk model (Levine, et al.,1987). The CEAREX data indicated 
a wavefield less energetic than lower latitude fields, but 
still significantly more energetic than the AIWEX wave field. 
The CEAREX data, due to location, were subject to large 
diurnal tidal motions which pumped energy into the internal 
wave field. In addition, under ice topography (ice keels) and 


local sea floor geography ( 1.e. Yermak plateau) affected the 
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internal wave field (Stanton). In contrast, the AIWEX data 
were taken in the Canada Basin of the Arctic Ocean (Padman, et 
al.,1990) and experienced far more quiescent conditions. 

Arctic internal wave fields possess different 
characteristics from lower latitude fields. Some suggested 
causes in the Arctic are less forcing due to wind, a damping 
effect caused by the ice cover, lack of surface waves, and 
weak large scale circulation (Levine, et al.,1987). Choice of 
frequency dependence in the internal wave field model 
critically depends on the location that is to be modeled. An 
arctic model will use a 1/f dependence for the velocity power 
spectrum, whereas a lower latitude model will use Ue The 
CEAREX location represents a special situation where large 
tidal forcing causes the internal wave field to possess 
characteristics similar to lower latitude fields. 

Garrett and Munk (1971) show a collection of data on 


displacement power spectra versus wavenumber obtained from 


towed measurements. The data follow a 1/k? pattern over a 
full four decades. The largest wavenumber corresponds to a 
wavelength of 10,000 nm. Assuming a simple plane wave 


description, the velocity spectrum will possess the same 1/k° 


dependence 


d(displacement) _ ; C eilkx-wt) 


displacement = —G_eitkret) | velocity = eT: 10 —— 


A simple empirical model for the k dependence could resemble 
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the omega dependence model with a power limiting ke factor, 
However, to facilitate later mathematical computations, use 


instead 


Ce 
Iv|? = on kan < k < @. (57) 
The total energy in the spectrum is now limited by the cutoff 
at k,,, instead of having the spectrum level out at low k. 
Assume this model applies to both arctic and lower latitude 
internal wave fields. 


Combining both forms, a model for the velocity spectra for 


the Arctic (i.e. AIWEX) 1s 


Ge e- 
alc ee es ee ae Sy enn k< ~. (58) 


Jwrt+w? k? Or eee ke 


The model for lower latitude and CEAREX velocity spectra 


becomes 


ee ee 
|v? = —————_. = ——__+___, Ie ge hie 59)) 
(Wo + w*) k? ANCE ar) ke 
Comparison with the AIWEX and CEAREX data determine the 
Constants of proportionality, C, and C,. “At 1 cph, the AIWEX 
spectrum is approximately 0.3 (cm/s)*/cph (0.11 (m/s)*/Hz), and 
all the CEAREX spectra have approximately the value 0.5 
(cm/s)*/cph (0.18 (m/s)*/Hz). Approximating the AIWEX model 


as proportional to 1/f at this frequency 


a7 


2 
oF ar ee oF 











Liers,: Serie err (60) 
Evaluating at 1 cph 
Ral Goel |P = 0. 3fen/s)” = aT’ (61) 
Then 
Ci = 2n(0.3)0Mem/s)? kk, = 1.9x20 & ee (62) 


Ss? 


Approximating the CEAREX model as proportional to 1/f* at this 


frequency 


|v.|? = ee ce (63) 
At? peo ee A fae Ie 


Evaluating at 1 cph 


Aem/s)? _ Ce 


Vell Gpay | —a0n5 IS alice 
Wc sl cph Ayla (ahora Ne. 


(64) 


Then 


2 2 m 
Co = 4n°(0.5) (cm/s)*-cph ann = 5 See 2 ean (—3? . (65) 
The formulas for the induced magnetic field involve the 
velocity amplitude and not the velocity spectral density. The 
conversion to a velocity amplitude was made using the 
definition of the power spectral density. The conversion for 


this case was simple because of the simple time dependence. 
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Given a function f(t), the power spectral density is the 
BouGier transform of tne autocorrelation function. of f(t). 
Let F(f) be the Fourier transform of f(t). The correlation 


theorem says (Brigham, 1988) 


” 66 
Beuroecorreiation of f(t) = { f(t) f(t+t) dt + FT pair ~ F(f) r( f° 


The power spectral density is thus F(f)F (f). Assuming the 


velocity has the following time dependence 


Vie Useeit) = Geet, (67) 


then the Fourier transform of the velocity is 


Ca ET Ce a ae Aloo No Ga gh) (68) 


The power spectral density becomes 


VAG ear) Vee ey = Oa ho) = Ge (69) 


Thus, for this simple time dependence and real coefficient C, 
the power spectral density is obtained simply by squaring the 
coefficient. 

Applying this result to the above model, the velocity 
becomes 


ae a, ae Uc) = eke a0? 


c 
Ne 4 Ww?) 1/4 k ly’, . PRY: k 
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This simple model in which the frequency dependence factor and 
the wavenumber dependence factors are multiplicative factors 
makes the following calculations much easier. However, as 
Garrett and Munk (1971) point out, actual observations 
indicate that a given wavenumber's contribution depends on the 
frequency. Following Garrett and Munk, the function u(w) was 
included here as a reminder of this additional frequency 
dependence 


Cy le vetkiw) = el) ge ee OM 


5 c 
(ie toa oe cease 


v,(k,w) = 


Equation [71] replaces the constant "C" in the induced 
magnetic field expression, Equation [50}. 

Examine two different methods to remove the k-dependence. 
For the first method, assume internal waves exist with a 
continuum of frequencies from zero to the buoyancy frequency, 
N, and assume they possess a continuum of wavenumbers from a 


minimum, k 1eorg slip y obo ete The rationale is that actual 


min! 
internal wave fields result from a complex superposition of 
waves generated by different sources propagating in a non- 
uniform ocean. The assumption of a continuum may circumvent 
a detailed description of this wave field. The second method 
utilizes the eigenfunction mode structure of the wave field 
required by Equations [17] and [19]. Application of a 


boundary condition on the vertical velocity at depth D yields 


a dispersion relation between w and k, for each mode. These 
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dispersion relations remove the k, dependence in B'. The 
power spectral density B'* is calculated for the lowest order 
mode n=1. The lowest order modes will contain a majority of 
the energy in an internal wave field. 

Try the continuum approach (first. TO» ‘simp hi fy 
Computations, assume that the depth D >> 1/k.. Then the last 
two terms in the induced magnetic field expression, which are 
weighted by exp(-k,D), can be neglected with respect to 1, 
simplifying the conversion to power spectral densities. The 
induced magnetic field has the simple e'*t time dependence 
referred to above. Squaring the amplitude gives’ the 


corresponding power spectral density 


ak 2 2 2 Z 
IBA? = [pBoc,]? 1 Ee) ee (72) 
Kx Niywotw? 
Z Z 2 2 
|Bi,|? = [p BoC]? — (w) we (N“+w*) (73) 


tee N4 (w+?) 


The formulas above are "one-dimensional" formulas in that 
they were derived assuming only one-dimensional wavevectors, 
k,. The more general wavefield will be two-dimensional with 


x 


wavevectors of magnitude 
k = [k2+k?. 78) 


The one-dimensional expressions for the induced magnetic field 
power spectrum must be converted to a two-dimensional form. 


Yaglom (1987) shows a way to accomplish this. Given a one- 


a4 


dimensional power spectral density function, f,(k,). Then the 
corresponding two-dimensional power spectral density function 
is 

r Ge LE (kee) dks 


k dk, 2 (75) 


Pai) See 
rs ki = ee 


The one dimensional induced magnetic field k dependence was 
‘We Appendix D contains a summary of the details in 


computing a two dimensional k dependence. Then 





mes ilk 
UR) = ri oe (76) 
Now, perform the integration over k (Keo < k < o) 
” 3 1 2 eb co 3 1 
a= eee dk = — eee eee a —————. 2 
J. aie 16 k4 Mun = 76 cn (77) 


The induced magnetic field two-dimensional power spectral 


density becomes 


z 2artve Zz 
Bh? = 1apeeen a ee ea eae (78) 
6 Knin N4 Vw o+ w? 
a Core 2 
|BoI? a [p .o BC]? a (Ww) (N*-0*)@ (79) 
4 4 2 2: 
Sa cee N* (W%+W*) 


Now apply the second method using the modal structure. 
Combining Equations [26],[{27] and [28] gives the vertical 


component of velocity 


Wee Cea ce Rem caine (for ky=0). (80) 


oe 


eu sew 6. — e721.) tien 


Ye ola oe) ome CEE (81) 


As a boundary condition, require that the vertical velocity 
equal zero at the lower boundary of the internal wave field at 


depth D. Then 


wien = C sin(6D) 6 =~ ° ms 0, (82) 
Zane 
sin(6D) = 0, => 6D = N*-0* k — 71T y jG a eee (83) 
G)* 
2r2 az2 
peat at | eer 2 _ Mee (84) 
= D Ze De? a7. 
en sath Soa at eB, 


Replace Kin the formula for B!' in Equatvon 9§[50] >with. this 


dispersion relation. Taking h=0 and mode number n=1 gives 





BS G) 
—2. (N?-w?) aes Neer ie (85) 


B= eosere 7¢2 
TN 


Converting to a power spectral density 





G) 
2 -T 
BS [p 0B] 2c? 2 (N2-w?)? ieee N?-w? ]2. (86) 
Toe 


Model the velocity spectra with only a frequency dependence 


Since a dispersion relation exists relating k, to o, 
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C? 
LA = Ivcl? = - 


(87) 
2 2 2 + Ww? 
Wi + w Oo * @ 


Replace C* in the B'* expression, Equation [86], with Equation 


[87] 





2 2 Sf | ae 
Bf? = ty o5)?@A_ 2 (newt)? [a+ e VO 92, (88) 


274 
w2+w? % N 





= “7 aes 
|B |? zs [p oB)?2 <2 — D* (N*-w*)? [1 + e N-w? 12, (89) 


w+ 0? 1™°N4 


C. UNDER ICE TURBULENCE MAGNETIC SIGNATURE 

1. Geometry 

The same geometry as for the internal wave calculation was 
used, except that now the y-coordinate was important (Fig. 5). 


The distance R was given by 


IR| = "Vee + oye Ge eee, with R= xf + Ve + ok. (90) 


For calculational purposes, turbulence was assumed to be 


Coniined to Ehe Teqiven Dope 
2. Turbulence magnetic signature 
The expressions for the earth's magnetic field and the 


induced current density are 
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field point B 


X 





O 
: Lt 
source point 
Natt AL dl eee Poel) NANCE eV |b) 


Z R = x + y+ (h+z) 


Rugure. 5: Geometry for Calculation of Turbulence Induced 
Magnetic Field 
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B=BK, and J=a(v~x B). (Se) 


Due to the vertical direction of the magnetic field in arctic 


regions, the current density has only x and y components 


J = o(v,Bl - v,Bj). (92) 


Make a change of variables for the vertical direction 


Z2+h=Z2;AZ—02; Wange DCzGme= Deen ae ae (93) 


Then 


S t 
Jx R= oBl-v,zi - v zi + (v,y+v,x) kK). (94) 


The Biot-Savart law then gives 


So a ee Bo _-v,zt - v,zt + (vyy+v,x) K 
—— d’r= 6B 
4nxn R3 4n = 


(x2 + y2 + 2?) 2 


dB! = 








dxdydz. (95) 


Integration over x,y,Z gives the components of the total 


induced magnetic field. 


i(k x-ot) z 


B = = = 
Bas ff (doen. va calc, fa eee ee EEE Gia? Gryhaan 


dxdydz. (9) 





The y-component equation is similar. The z-component is 


ele | [fates 





xX=ue ine | Kx ik,y Vie e ive 
a= f’ Ae Slee JXCk eee eg re 


~Dosh i [kz=k? otro? (x?+y2+22) 2/2 


Performing the integrations [see Appendix B]) yields 


a6 


-jwt 
i. alee Oe = eq let Be [oe *ePe - o *x?) 
x 


2k ko+kx (wo >+ ow? 


(98) 


= .0oByle.?* : ; 
pi = — Po PY git ety panko _ g7hy?y | (99) 


2k ki+k; yw+ 0? 


<< 


ip ,oBv,e77°° olen ene. coe een ene) ) , (100) 


2,/w%,+ 0? fonieeeice kywWko+ky 


BL = 


Approximating the geometry for an on-the-ice measurement 


6 Maem Or D, = 0; D =o, (101) 
-p oBve 7° 
Bee = ee —— (102) 
2 W ot Ww? Ky Ko ey) 
fod TBOB YE 7S" 


B, = (103) 


2/w4+w2 Kal ki+k? ie esa ee 


Letting D approach infinity greatly simplifies the following 
calculations. This procedure is justified due to the expected 


magnitude of the parameters 





k, = 22 ~-2%_;  pe100-1000m; =e < 0.002 < 1. (104) 
A 100m 
These equations still have a k-dependence. To get the 


components as functions of only w and t, integrate over the k- 


dependence 


a4 


-p oBv.e7°° pk =0 il 
Bi = ef —_———. dk, , (105) 
2 Oot O > ie iene) etc 
-1 oBv.e 7% k2+k2., + k 
_ Hot BV oe ane may) (piccecem eon (106) 


2 w*+w? IS in 


where k.. is the smallest wavenumber (corresponding to the 
longest wavelength) present in the x-direction. A limit on 
wavelength is realistic for the ocean, and required to keep 
the integral finite. Recall that k, keeps the total power 
finite as k approached zero. k, thus represents the poiniaam 
the power spectrum at which the characteristic k* behavior 


begins. k, is thus a reasonable choice fork... Set k- Sequam 


n 


to k.- Then 


-L oBv.e ie 
Bi = Pi Sone OW aha an ive th) | (107) 


20 +0? o 
Similarly for the y and z components 
Bl = -p ,oBv,e 77° LES — Ble ip,oBv,e7%** in YE +4) (108) 


20 %+o? ° yo oto? ° 


For the z-component, integration was done only over the k 
components which gave rise to the respective term via the 


Biot-Savart formula. 
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The total induced magnetic field is now 


oBy.e77°* 
Bo O Anty2 +1) Jé . (109) 


270 o+ 0? ° 


B/ a Se Bie + ish - IS ioe 


The horizontal component of the induced magnetic field is 


Pe ee ee Nike 
Bi, = Be +B a ae 2 SS SS Sn tvesh) (110) 


V2 YW ot W? o 


Again, the power spectrum was obtained by Fourier 
transform of the autocorrelation of B', which is equivalent to 


squaring the amplitude here 


Bl? = 3 (pH oBv,)* (In(y2 +1) )* (111) 
Die). 


(ie By.)° Cink 241) )- 
Bi’, = : “ot ee; a (112) 


D. MEASURED MAGNETIC SPECTRA 

The Arctic Internal Wave Experiment (AIWEX) was conducted 
Moo OlstLhe arctic 1Cce at 74°N~ b45°W. ~ in April of that 
year, P. Czipott and W. Podney used two induction coil 
magnetometers to measure the North-South horizontal component 
of the Earth's magnetic field. Fig. 6 shows the results from 
three runs. The spectra flattened out at low frequencies 
because the electronic gain declines at low frequencies. The 
region in the 10° -10% Hz decade was relatively unaffected. 


The spectra here showed a characteristic 1/f* dependence. 
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They also showed an approximately hundred-fold magnitude 
variation between the three runs. Thus, large magnitude 
changes can occur over a temporal scale on the order of days. 
Note that the bottom spectrum basically coincides with the 
dashed line spectrum which corresponds to measurements made on 


land at Poker Flat, Alaska. 
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Power spectra for magnetometer output from Runs 100 
(center), 103 (top), and 107 (bottom) at the ice camp, 
and Run 124 (dashed curve) from Poker Plat. The scale 
is calibrated by the magnetometer gain at 0.1 &z: 
Pigure 3 illustrates the dropoff in gain at lower 


frequency. Error bar shows the 90% confidence 
interval. 


Figure 6: AIWEX Power Spectra for magnetometer (Czipott, 
Podney, 1985) 
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IV. RESULTS 


A. EVALUATION OF INTERNAL WAVE MAGNETIC FIELD MODEL 
Numerical evaluation of the formulas from the model used 


the following values/parameters 





of Bed 7 S, 7 , 
LL, = 4x10 a Oo = ao eee B= 25,000 nT; (is) 
9P -3 33x70 kg. N=1.05x10-2224, (114) 
OZ m4 sec 
iene = LO mM; = Keet= 6 28x 10- ne (115) 


Table II lists sample numerical values using the continuum o- 
and k-dependence method. Table III lists sample numerical 
values using the modal method. For the continuum method, Fig. 
7 shows the AIWEX model spectrum, the CEAREX model spectrun, 
and an extrapolated spectrum for the AIWEX magnetometer data. 
Extrapolation was made from the 10°' - 10% Hz range of Run 
100. Fig. 8 shows the same spectra using the modal method for 
the AIWEX and CEAREX models. The vertical extent, D, of the 


internal wave field was 300 m (D'Asaro, Morehead, 1991). 
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Figure 7: AIWEX Model Spectrum, CEAREX Model Spectrum, and 
Extrapolated Spectrum for AIWEX Magnetometer Data (Continuum 
Method) 
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Figure 8: AIWEX Model Spectrum, CEAREX Model Spectrum, and 
Extrapolated Spectrum for AIWEX Magnetometer Data (Mode 
Structure Method) 
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Table II: Internal Wave Model Numerical Values 
(Continuum Method) 





frequency (Hz B' AIWEX (nT¢/Hz) B' CEAREX (nT°/Hz) 





10° 1.3x10° 6.0x210° 
10"? boo" 6.0x10° 
10° B21 0" 3.8x10° 
1.6x10° 1s Os 4.9x10° 
Table III: Internal Wave Model Numerical Values 


(Mode Structure Method) 





Frequency (Hz B' AIWEX (nT*/Hz) B' CEAREX (nT‘*/Hz) 


10° DSO. sense ills 
Ho: DE OF 1.3x10° 
10° dele Gs 5.2x10° 
1.6x10° 1.2x10° Jeo <0; 


The AIWEX magnetometer measurements were not compensated 
for ionospherically generated magnetic fields. To separate 
the internal wave induced magnetic fields from the ionospheric 
fields, Czipott and Podney placed a tiltmeter on the ice as an 
independent measure of internal wave activity at the AIWEX 
Site. They also obtained magnetograms of ionospheric activity 
during the experiment from several measurement stations in 
Alaska. Following Czipott and Podney, the signals from these 


three sources are compared in Table IV. 


Table IV: Signal Characterizations 
Bc Me ee em ee, a eee 


RUN # IONOS PHERE TILTMETER MAGNETOMETER 
(MAGNETOGRAM) 

100 "quiet" "high" "active" 

103 "active" wactive. "high 

105 "active" Hoque! "quiet" 

107 "oulet" "quiet" "quiet 


As a first order analysis, assume the lonosphere had only two 
levels of activity, "quiet" and "active", measured by the 
magnetograms. Assume the under ice seawater motion also had 
only two levels of activity, “high"/"active” and “quem 
measured by the tiltmeter. Then the above four runs show the 
resulting magnetometer signal for the four’ possible 
combinations of lonospheric activity and seawater motion. 
Runs 100 and 105 indicate that the lonosphere and the seawater 
motion each produced a measurable magnetometer signal. 
Together, runs 100, 103, and 105 indicate that the magnitudes 
of the magnetometer and ionospheric signals were of the same 
order. Run 107 affirms that the magnetometer readings result 
from actual ionospheric and seawater motion signals. 
Extrapolation of the AIWEX magnetometer data with a 1/f? 
frequency dependence to the 10> = (lok ez region gave power 
spectrum values between 10° and 10° nT*/Hz. For the continuum 
method, the model values varied one to two orders of magnitude 
above and below these extrapolated values. The 17 


dependence made these continuum models sensitive to the value 
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ef k. chosen. Choosing A... = 1000 m instead of 10,000 m 
reduced the values by an order of magnitude. For the modal 
method retaining only the lowest order mode, the AIWEX model 
predicted values one to two orders of magnitude higher than 
the extrapolated data, while the CEAREX model predicted values 
about two orders of magnitude above the extrapolated data. 
The mode structure model values depended on vertical extent of 
Gherwinternal wave field as D*. Taking D°="100 m instead of 
300 m reduced the spectra by one order of magnitude and 
brought them within one order of magnitude of the extrapolated 
data. The vertical velocity profile of Equation [{81] has a 
Simple sin(éz) dependence. Refining this profile could 
significantly reduce the energy in each mode since B!? is 
proportional to (velocity)®*. Overall, the models reasonably 
predicted the magnitude of the seawater induced magnetic 
fields. 

The frequency dependencies of the continuum models varied 
significantly from the 1/f* dependence of the extrapolated 
data. The AIWEX model predicted a linear increase with 
frequency and then a rapid falloff as the buoyancy frequency, 
N, 1S approached. The CEAREX model predicted a parabolic 
falloff with frequency as (N*-w*). Unfortunately, the gain 
falloff of the AIWEX data over internal wave frequencies may 
have masked any dropoff near the buoyancy frequency. These 
continuum models predicted too much energy at the higher 


frequencies. Apparently the assumption of frequency and 
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wavenumber continua lead to erroneous’ predictions. The 
specific relation between w and k, denoted by p*(w) in the 
models, must be determined. If one assumed that u(w) goes as 
1/(number of wavenumbers), and that number of wavenumbers was 
proportional to frequency (see, for example, Garrett and Munk, 
1972, for a similarly modeled dependence), then the models 
would more closely predict the w* dependence of the data. One 
could compare the internal wave model in the 10°! - 10% Hz 
range because internal wave oscillations are limited to 
frequencies below the buoyancy frequency. 

The mode structure method gave models with a much more 
realistic frequency dependence. The AIWEX model declined as 
1/f before going to zero near the buoyancy frequency. This 
dependence paralleled the observed AIWEX velocity spectra 
frequency dependence. The CEAREX model fell off as 1/f*, just 
as the CEAREX velocity spectra did. These parallel frequency 
dependences between induced magnetic field and generating 
velocities satisfied physical intuition. This close 
correspondence between model and data implied that a proper 
description of internal wave induced magnetic fields requires 


utilizing the internal wave modal structure. 


B. EVALUATION OF TURBULENCE MAGNETIC FIELD MODEL 
Numerical evaluation of the formulas use the following 


values/parameters 
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po = 4pixi07™; o=4.0 2; B= 25,000nT; (116) 
A? m 
3 a m? = -4 
Vo = 5.7xX10°4 k, ——; Ge = et eC i (117) 
Hz-s‘4 


The Wee value was determined the same way as Ce=fer the 
internal wave model, assuming that the CEAREX data 1/f? 
dependence extrapolates to the 10°'-107 Hz range. Sample 
numerical values using Equation 102 are given in Table V. 
pide? sshows a plot cof the B,* power spectral density (nT*/Hz) 
versus frequency for Aes = 100 m and also the extrapolated 
AIWEX magnetometer spectrum. 


Table V: Turbulence Model Numerical Values 


FREQ (Hz) B'.?(nT*/Hz) (A, =100 m) 
10° 4.4x10° 
10°¢ 4.4x10° 
5x10° janice be 
TO 4.4x10°'° 


mies medel. spectra, exhibited.a characteristic. 1/u° 
dependence, just as the AIWEX data show in the 10°' - 10° Hz 
range. However, the model values were seven orders of 
Magnitude less than the AIWEX data. At 10° Hz the model 
predicted magnitude was approaching the limits of the most 
modern field detectors at 10° (nT*/Hz) . The smallness of the 
predicted signal derived from the small velocities used and 
the maximum characteristic wavelength of 100 m. The small 


velocities on the order of 10°' - 10% (cm*/s*-Hz) were 


experimentally measured during AIWEX (McPhee, 1989) and their 
power spectrum is shown in Fig. 10. The assumption of 
isotropic turbulence in the model restricted wavelength values 
to less than about 100 m because the maximum observed depth of 
the upper boundary layer is about 100 m (Stanton) (note: the 
"average" depth is about 30-40 m). Thus, while turbulent 
motion in the upper boundary layer of the ocean could produce 
magnetic signals, the magnitude remained small because the 
volume of water generating the signal was limited. This 
implied that the AIWEX magnetometer measurements in the 107! - 


10°°> band were due to ionospheric signal only. 
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Figure 9: Turbulence Model and Extrapolated AIWEX 
Magnetometer Spectra 
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Power spectra of the vertical] velocity ume series at each level. The 
95 percent confidence interval is shown at lower left. The spectra are 


plotted on log-log axes, with successive spectra displaced one decade down 
for clarity. 


Figure 10: AIWEX Vertical Velocity Power Spectra 
(McPhee, 1989) 
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V. CONCLUSIONS 

The magnetic fields induced by the motions of seawater in 
the earth's magnetic field were modeled analytically. Such 
analyses required mathematically modeling the velocities of 
the internal wave field and internal turbulence field, and 
then applying Ampere's law via the Biot-Savart formulation. 
For internal waves, the values predicted by the model for the 
induced magnetic field were roughly equal in order of 
magnitude to extrapolated magnetometer data. Treating 
internal wave frequency and wavenumber as continua yielded 
unphysical frequency dependences. Using an eigenmode 
structure for internal waves yielded more realistic frequency 
dependences. Thus, proper analysis of internal wave magnetic 
Signatures must incorporate the internal wave modal structure. 
Additional magnetometer measurements should certify the 1/f° 
frequency dependence and determine the behavior near the local 
buoyancy frequency. The turbulence model successfully 
coincided with magnetometer measurements in frequency 
dependence. The predicted frequency dependence matched the 
yee dependence of the data. However, the predicted 
magnitudes were several orders of magnitude below the AIWEX 
magnetometer measurements and approached the limits of the 
most sensitive magnetometers. This implied that the 


relatively small scale of turbulence motions (A™* = 100 m) 
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resulted in a small induced magnetic field signature. The 
AIWEX data in the 10! - 10% Hz range were interpreted as 


deriving from ionospherically generated signals . 
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VI. RECOMMENDATIONS 


Additional measurements of seawater motion induced 
magnetic fields are needed. Such measurements should be made 
with hagh ly sensitive magnetometers, fOr example, 
Superconducting quantum interference devices (SQUIDs) . 
Discrimination against lionospherically induced magnetic fields 
must be done by using the large spatial coherence length of 
ionospheric fields to identify and remove their signal. 

The models can be extended to predict the signal seen by 
a gradiometer. This involves expressing the induced magnetic 
field at the field point in terms of field point position and 
taking the gradient. Gradient measurements of internal wave 
induced magnetic fields are more numerous than magnetometer 
measurements (for example, Czipott and Podney, 1985; Podney, 
1975). 

Both the internal wave and the turbulence models for 
induced magnetic field utilized the approximately vertical 
orientation of the earth's magnetic field in the Arctic to 
simplify the calculations. The models can be extended to 
apply the Biot-Savart calculation method to the non-vertical 
magnetic field at other parts of the world's oceans. 

The internal wave model assumed a constant buoyancy 
frequency, N, for the entire wave field in the horizontal and 


vertical dimensions. In the actual ocean, the buoyancy 


5S 


frequency changes from location to location and varies 
significantly with depth. A more refined model must include 
such variations and should predict a less defined cutoff 
frequency due to the spread in buoyancy frequencies. 

The modal version of the internal wave model should be 
refined by including the contributions of higher order modes. 
The resulting correction should not be too large because 
energy content rapidly declines for higher order modes due to 
dissipation. The velocity profiles should also be refined. 
As an example, require velocity continuity at the boundaries 
of the internal wave field and allow exponential decay outside 
these boundaries, analogous to going from an infinite square 
well solution to a finite square well solution. The amplitude 


Watchineweme welleis seduced 
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VII. APPENDICES 


A. CALCULATION OF INTERNAL WAVE INDUCED MAGNETIC FIELD 
For a region with constant Brunt-Vaisala frequency, N, the 


vertical component of the seawater velocity is given by 


AO ra ara are an ial) va Sees (118) 
WZ) waen ea gnGie fen Gren ae. (119) 
N2-w? 
where: 5 =k, ,| ———., k,=0. (120) 
G)° yi 


Require aS a boundary condition 


Mae oO ee) 0, ee Ow SHO) aCe tr C |\ae a2 Geese] 
Then the vertical velocity component becomes 
Wiz) = Cer = [els =e 207) (122) 


The derivative with respect to z is 


SS) ea Clete) me een toe. = =. (123) 


With the simplified geometry, the y-component of velocity is 


zero. The x-component of velocity becomes 


a7 


BUSING) Beis) OE is ~ kW) (==) BESS © ie) (124) 
Yh, 
= = [e2824 e@- 262] NS) . (125) 


x 


These expressions can now be used to calculate the induced 


magnetic field 





/ _ Ho, di 
By = iL " rr ee. (126) 





= +00 =D 0 wB = uB 
ee (a [i Bay) SEEPS al) spas az (127) 
x=-0 JZ=D, Za R? 
(e282. @-162) cosycosd + fas (@i'2, gibt) (128) 


sil By oc Wes tree? [h+z] dzdx 


© /z=D, 2n (h+z)?2 + x? 


The x-integration is done via contour integration in the upper 


half plane 
Seek ik ik ik 129 
f e x* 2 f e xX* ee Z f e x* ( ) 
x=-= = (h+z)?2 + 22 unp (x+1(h+z)) (x-1(h+Z)) bndry (h+z)*+z? 


ik, (i(h+z)) Spe (130) 
e 0 = T e k,(h z) ; 


a Sere grr ary rs ~ Chez) 


Then 
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ar ae 131 
Bin EOE gree tae f(g tril t 9 “He l)2) cosycoss : Ssi¢ (a (RE, 9 (ky 182) 7 ( ) 
De a 


Performing the z-integrations and noting limits of integration 





yields 
po. HB BoC “kh -iwt 1 (-k,+16)z 1 (-k,-ib)z (22) 
B, = mare ce f ome + ( mate ye }cosycos® + 
dsing | 1 \ @ (kat id) z_ 1 ye (ker i8)z7 te (133) 


Ps “k,+16 k,+16 


x 


The z-component of the induced magnetic field is 


calculated similarly 

















) _ pxtem peed *< 7D ee o (WB, ~ uB,) 

Be ff, gear * * [oh s2 me x ae (134) 
-fpeteees pees (@162_ @- 182) cosycosd + Being tellece-t) | et Frx-°2) aedx(135) 

ay (h+z)? tuixe 
2 eee ace (e#?-e°19) cosycoso+ SS1nP (1824 @-182) ] ime hethee) gz (236) 

ae Le oe ktig-iwt (-k,+18)z 1 (-k,-18)z 
B, ee =m ss + ( fn © Jeosycos@ +(137) 
, Ssind [( ye (kxed8)z_¢_ 1 ye (Rr ibz) 3 ee (138) 

ke k a8 k,+18 re 


Note that the x and z components have the same magnitude and 
differ in phase by 7/2. These last equations can be 
Simplified by multiplying the numerator and denominator by the 
complex conjugate, and using Euler's identity for complex 
exponentials, resulting in 


HP BoC ( uh =Ky (Nez) (139) 


+ 
<a 


iBl=Bl= ; ; 
+ 
xXx 


de eft? (-k sindz - dcos6z)2icosycosd + 
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+ (-k,cos8z + bsingz) 2°51ne ) [s.- (140) 


So 


B. CALCULATION OF TURBULENCE INDUCED MAGNETIC FIELD 
The integrations for the components of the turbulence 


induced magnetic field are shown below. 








HW oB y=@ ¢z=D+h -v,e tee e ik 
By = —__—_._.——- dxdydz, 
x 4 ey ial a +h Yk2+k foe +0? Gtayta gay iP a4 Zz (141) 
B px=@ pz=Deh =e yo ») dy 
[haat +h 2 2 im (x*+y7+27) 3/2 dydz, (142) 
: cet hee 
=o -iwt a eo ikxx 
p Bohs i1an 2, [rrr 2{% A ae (stm 
Ati tie OT O> aac ee xe ZG xe ep 
-2y ,oBv.e 7% z=Dth 2 @ ke 
i — ad “f = Se (144) 
xs-Wz=D+h (x*+7%) 


any ko+ké yw +w)* 


Perform x-integration in upper half of complex plane 


-2 OoBV =D+h _ 
Bi wai loader Tt is ex? QZ, (145) 
anykitke = “7 Porh 
oBV a = ~ 
2 ok 26 gk. eh eo eee (146) 


2k ylki+k? 


The integration for the y-component is performed in the exact 
same manner, with x and y interchanged. 
The z-component integration is calculated similarly 


ye 
ae we Aes +h etn ee +22) 3/2 er tics 


pt state 


B, 





1,(147) 
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-iet 
Leh y 2 ( 
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4nywi+w? 
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=< x etk“dxdz 
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ad x00 kisk? ((xt+y24z2) 3/2 


. a y ee" dydz 


nea 2ye*dxdydz 


1,(148) 


1,(149) 
a hgseo ce (y*+z22) 
y et aydz (150) 


sy k2+k? (y+iz) (y-iz) 


DRY. 
ky ikot+ky 


re 


Bee, 
Wutem a ea, 
yk pee 
) Seem een) (12) 
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INTERNAL WAVE VELOCITY POWER SPECTRAL DENSITY MEASUREMENTS 


from the Coordinated Eastern Arctic Experiment 


(CEAREX), March-April 1989. 
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Rotary Spectra 250m S4 
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D. 1-D TO 2-D POWER SPECTRUM CONVERSION CALCULATION 
The one-dimensional power spectrum and its derivative are 


Cees) 
Ek.) = ae a = coe (153) 
ak 


1 


The two-dimensional power spectrum is given by 


1 pe af, (k,) dk, 
(2 (ee is S 
A) tm Jk dk, pomine (154) 
1 
oe + sag ak, 
ee eee (155) 


i) oe 


Make a change of variable 





mee Pee <, dk, = -—<dn, (156) 
(157) 


ft 
Moke lo ftom? 


From CRC Standard Math Tables (1987), given a trinomial inn, 


M= a+ bm + cm*, then 








ges iW2med) ae (n-1)a fm"? 
an-= ni /M ee 
i IM ee ty lone i IM 2) 


Here, a=1, b=0, c=-1. Applying this integral twice yields 


Ce oS 
Pees oe eat ye oman 4 Ssintm lf, (159) 
tk? 4 8 
: ae 
fo.) = i ee (160) 
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